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We study spontaneous radiative decay of translational levels of an atom in the vicinity of a 
semi-infinite dielectric. We systematically derive the microscopic dynamical equations for the spon- 
taneous decay process. We calculate analytically and numerically the radiative linewidths and the 
spontaneous transition rates for the translational levels. The roles of the interference between the 
emitted and reflected fields and of the transmission into the evanescent modes are clearly identi- 
fied. Our numerical calculations for the silica-cesium interaction show that the radiative linewidths 
of the bound excited levels with large enough but not too large vibrational quantum numbers are 
moderately enhanced by the emission into the evanescent modes and those for the deep bound levels 
are substantially reduced by the surface- induced red shift of the transition frequency. 
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I. INTRODUCTION 

Due to recent progress in atom optics and nanotech- 
nology, the study of individual neutral atoms in the 
vicinities of material surfaces has gained renewed inter- 
est 0, 0j S 0, @- A method for microscopic trapping 
and guiding individual atoms along a nanofiber has been 
proposed 1] . An optical technique for loading atoms into 
quantum adsorption states of a dielectric surface has been 
suggested 0, 01 ■ The possibility to control and manipu- 
late individual atoms near surfaces can find applications 
for quantum information [H, 0, @ an( A atom chips [jj, • 
Cold atoms can be used as a probe that is very sensitive 
to the surface-induced perturbations Such state-of- 
the-art applications require the development of a system- 
atic microscopic theory for spontaneous radiative decay 
of atoms near surfaces. 

It is well known that the spontaneous emission rate 
of an atom is modified by the presence of an inter- 
face 0, 0, HH flEl fill . A quantum- mechanical linear- 
response formalism has been developed for an atom close 
to an arbitrary interface LL-1 ■ An alternative ap- 

proach based on the mode expansion has been used for 
an atom near a perfect conductor Jl5| or a dielectric ^(| . 
In the previous treatments 0> ll3lll4L , the effect 

of the center-of-mass motion of the atom on the sponta- 
neous emission process was neglected. In this condition, 
the effects of the surface on the spontaneous radiative 
decay of the atom manifest simply as a quantum electro- 
dynamic enhancement due to the mode modification and 
a frequency shift due to the surface-atom van der Waals 
interaction, with no structural changes in the dynamics 
of the process 0, 0] . 

Lately translational levels of an atom in a surface- 
induced potential have been studied 0, Q . An experi- 
mental observation of the excitation spectrum of cesium 
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atoms in quantum adsorption states of a nanofiber sur- 
face has been reported [ljj ■ In the theoretical treatments 
of Refs. t ne quantum electrodynamic enhancement 

of spontaneous emission due to the mode modification 
was completely neglected, and a set of equations with 
phenomenologically added decay terms was used for de- 
scribing the dynamics of the loading process. A more 
rigorous theory requires a deeper knowledge of the struc- 
ture of the decay equations and the magnitudes of the 
decay coefficients for the translational levels of the atom. 

In this paper, we study spontaneous radiative decay of 
translational levels of an atom in the vicinity of a semi- 
infinite dielectric. Our study is general in the sense that 
it incorporates the quantum nature of the radiation field 
in the presence of the dielectric as well as the quantum 
center-of-mass motion of the atom. We systematically 
derive the microscopic dynamical equations for the spon- 
taneous decay process, which reflect the complexity of 
the translational levels of the atom in the potential. We 
focus our attention on the radiative linewidths and the 
spontaneous transition rates. Our results obtained by the 
quantum mode-expansion approach allow us to identify 
the roles of the interference between the emitted and re- 
flected fields and of the transmission into the evanescent 
modes. 

The paper is organized as follows. In Sec. [H] we de- 
scribe the model. In Sec. lIIII we derive the basic dynami- 
cal equations for the spontaneous radiative decay process. 
In Sec. lIVI we study analytically the decay rates. In Sec. 
[V]we present the results of numerical calculations. Our 
conclusions are given in Sec. IVII 



II. MODEL 

We consider a space with one interface. The half-space 
x < is occupied by a nondispersive nonabsorbing di- 
electric medium (medium 1). The half-space x > is 
occupied by vacuum (medium 2). We examine an atom, 
with an upper internal level e and a lower internal level 
<7, moving in the empty half-space x > 0. 
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A. Quantum translational states of the atom 

We describe the quantum translational (center-of- 
mass) states of the atom by using the formalism of Rcf. 
@. The potential energy of the surface-atom interaction 
is a combination of a long-range van der Waals attrac- 
tion — C3/X 3 and a short-range repulsion 0. Here, C3 
is the van der Waals coefficient. We approximate the 
short-range repulsion by an exponential function Ae~ ax , 
where the parameters A and a determine the height and 
range, respectively, of the repulsion. The combined po- 
tential depends on the internal state of the atom (see Fig. 
Uj l. and is presented in the form 



where j = e or g labels the internal state of the atom. 



(1) 



Here, p and m are the momentum and mass of the atom, 
respectively, and u)j is the frequency of the internal level 
j. We introduce the notation \<p a ) and \(pb) for the eigen- 
states of the center-of-mass motion of the atom in the 
potentials V e {x) and V g (x), respectively. The wavefunc- 
tions if a (x) and <pb(x) are determined by the stationary 
Schrodinger equations 



2m dx 2 



2m dx 2 



V 9 (X) 



<fa(x) = S a ip a (x), 

ipb(x) = £bipb(x), (3) 



respectively. The eigenvalues £ a and £b characterize the 
energies of the translational levels of the excited and 
ground states, respectively. Without the loss of gener- 
ality, we assume that the center-of-mass eigenfunctions 
ip a and ifb are real functions. 
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FIG. 1: Surface- induced potentials for ground- and excited- 
state atoms. The parameters are Cz g = 1.56 kHz |im 3 , Cz £ = 
3.09 kHz ^m 3 , A g = 1.6 x 10 18 Hz, A e = 3.17 x 10 18 Hz, and 
ct g — a e = 53 nm -1 . They are taken for the case of a D2-hne 
cesium atom in the vicinity of a semi-infinite silica medium. 



The potential parameters C33, Aj, and ay depend on 
the dielectric and the atom. In numerical calculations, we 
use the parameters of fused silica, for the dielectric, and 
the parameters of atomic cesium with the D2 line, for the 
two- level atomic model. In Appendix ^ the parameters 
of the ground- and excited-state potentials for the silica- 
cesium interaction are estimated to be C^ g — 1.56 kHz 
^m 3 , C 3e = 3.09 kHz ^m 3 , A g = 1.6 x 10 18 Hz, A e = 

1 The potential 



3.17 x 10 18 Hz, and a g ~ a e — 53 nin 



depths are D g = 159.6 THz and D e = 316 THz. The 
minimum positions for both potentials are assumed to 
be the same, equal to x m — 0.19 nm. 

The Hamiltonian of the atom in the surface-induced 
potential is given by 



FIG. 2: Energies and wave functions of the center-of-mass 
motion of the ground- and excited-state atoms in the surface- 
induced potentials. The parameters of the potentials are as 
in Fig. Q The mass of atomic cesium m = 132.9 a.u. is used. 
We plot only the bound levels of the ground and excited states 
with energies in the range from —1 GHz to —5 MHz and the 
free ground-state level with energy of about 4.25 MHz. 
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FIG. 3: Eigenvalues for the bound levels of the ground and 
excited states of the atom in the range from —1 GHz to —20 
kHz. The parameters used are as for Fig. [5] 



FIG. 4: Same as Fig. [5]but for the twenty lowest levels of the 
excited state and the forty lowest levels of the ground state. 
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(2) We introduce the combined eigenstates |o) = |e) <g> \tp a ) 

and \b) — \g) ® \<fb)j which are formed from the internal 
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FIG. 5: Same as Fig. [3]but for the twenty lowest levels of the 
excited state and the forty lowest levels of the ground state. 



and translational eigenstates. The corresponding ener- 
gies are Tiu> a = Tioj e + £ a and hujb — huj g + £b- Then, we 
can represent the Hamiltonian (J2J) in the diagonal form 
Ha = J2 a hu a \a){a\ + J^b ^b\b)(b\. We emphasize that 
the summations over a and b include both the discrete 
(£ a ,b < 0) and continuous (£ a ,b > 0) spectra. The levels 
a with £ a < and the levels b with £ & < are called the 
bound (or vibrational) levels of the excited and ground 
states, respectively. In such a state, the atom is bound to 
the surface. It is vibrating, or more exactly, moving back 
and forth between the walls formed by the van der Waals 
potential and the repulsion potential. The levels a with 
£ a > and the levels b with £b > are called the free 
(or continuum) levels of the excited and ground states, 
respectively. The center-of-mass wave functions of the 
bound levels are normalized to unit. The center-of-mass 
wave functions of the free levels are normalized per unit 
energy. For bound states, the center-of-mass wave func- 
tions ip a and ipi, can be labeled by the quantum numbers 
v a and fb, respectively. For free levels, the conventional 
sums over a and b must be replaced by the integrals over 
£ a and £b, respectively. 

The ccntcr-of-mass wave functions with energies in the 
range from —1 GHz to —5 MHz are plotted in Fig. The 
center-of-mass eigenvalues in the range from —1 GHz to 
—20 kHz are displayed in Fig. The mass of atomic 
cesium m = 132.9 a.u. is used in our calculations. Note 
that the energy levels in the above range are very high 
compared to the depths D e = 316 THz and D g = 159.6 
THz of the potentials V e and V g , respectively. Therefore, 
the corresponding wave functions are mainly determined 
by the long-range van der Waals interaction |2( ■ We find 
that the maximum quantum numbers for the excited- 
and ground-state potentials are 437 and 311, respectively. 
When the vibrational quantum number v is small, that 
is, when the translational energy level £ v is deep, the 
wave function ip v depends strongly on the short-range 



repulsion (see Fig. 0}. For such bound states, the center- 
of-mass wave functions are well confined in the proximity 
of the surface, and the surface-induced frequency shifts 
of the atomic transitions are substantial compared to the 



optical frequency ujo = ui e — 
D 2 line (see Figs. g|and[5J|. 



= 352 THz of the cesium 



B. Quantum radiation field and its interaction with 
the atom 

We use the formalism of Ref. 0] to describe the quan- 
tum radiation field in the space with one interface. We 
label the modes of the field by the index /.i = (u>K.qj), 
where tu is the mode frequency, K = (0,k y ,k z ) is the 
wave- vector component in the surface plane, q = s,p 
is the mode polarization, and j = 1,2 stands for the 
medium of the input. For each mode fi = (wKqj), the 
condition K < krij must be satisfied. Here, k — u>/c is 
the wave number in free space, rt\ > 1 is the refractive 
index of the dielectric, and n 2 = 1 is the refractive index 
of the vacuum. The mode functions are given, for x < 0, 
by 



U^Ksl^) 

U wKp i(a;) 
U wKs 2(a;) 
U W K P 2(aO 



and, for x > 0, by 
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In Eqs. {U and ©, the quantity j3j = (k 2 



K 2 ) 



(5) 



2U/2 



with Re (3j > and Im f3j > 0, is the magnitude of the 
x component of the wave vector in the medium j. The 
quantities rfj = (Pi-/3j)/(Pi+Pj) andtf^- = 2/3 i /(^ i +/3 j ) 
are the reflection and transmission Fresnel coefficients 
for a TE mode, while the quantities r v - = (Pin 2 — 
P 3 n 2 )/{p in 2 + P 3 n 2 ) and t£ = 2n l n J P l /(p z n 2 + P 3 n 2 ) 
are the reflection and transmission Fresnel coefficients 
for a TM mode. The vector s = [K x x] is the unit vec- 
tor for the electric field in a TE mode, while the vectors 
Pj'± — (KSc =F PjtK.)/knji are the unit vectors for the 
rightward- and leftward-propagating components of the 
electric field in a TM mode in the medium j' . 

Note that a light beam propagating from the dielectric 
to the interface may be totally reflected because n\ > 
112 = 1. This phenomenon occurs for the modes [i = 
(wTCqj) with j = 1 and k < K < kn\. For such a mode, 
the magnitude of the x component of the wave vector 
in medium 2 is P% = i\J K 2 — k 2 , an imaginary number. 
This mode does not propagate in the x direction in the 
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vacuum side of the interface but decays exponentially. 
Such a mode is an evanescent mode. 

The total quantized electric field is given by 0| 




'U^OKe-^'+H.c., (6) 



where a M is the photon operator for the mode p, 

and E P = E«p Jo°° J nj KdK ^ d(f> is the gen- 
eralized summation over the modes. Here, <fi is the 
azimuthal angle of the vector K in the yz plane. 
The commutation rule for the photon operators is 
[a M ,aL] = S(uj — u>')5(K ~ K.')8 qq >8jj> . When disper- 
sion in the region around the frequencies of interest is 
negligible, the mode functions satisfy the relation 
jT^dxn 2 {x)\Jl^ q .{x)\5 UJ ^ ql j l {x) = 2irc 2 (f3 J /uj)5{LJ - 
Lu')5 qq i5jji . Here, n(x) — n\ for x < 0, and n(x) = n 2 
for x > 0. Hence, we can show that the energy of the 
field is e / drn 2 (x)\E{r)\ 2 = £ M (ata M + o M oj,)/2. Here, 

[ dr = f°° dx f°° dy f°° dz is the integral over the 
whole space. 

We now present the Hamiltonian for the atom-field 
interaction. In the dipole and rotating-wave approxima- 
tions and in the interaction picture, the atom-field inter- 
action Hamiltonian is 

H iat = -ih G^aba^- 1 ^^ 1 + H.c. (7) 

fiab 

Here, a ab = l a )(^| describes the atomic transition from 
the level b to the level a, Lu ab = to a — cu b is the angular 
frequency of the transition, and 

G, iab = _— i_ e jK - R (tp^ • d e9 |^ b ) (8) 

is the coupling coefficient. In expression (jSJ, d eff is the 
dipole moment of the atom, and R = [y, z) is the projec- 
tion of the position vector r = (x, y, z) of the atom onto 
the interface plane. 



Meanwhile, the Heisenberg equation for the photon oper- 
ator a M is a M = Y,ab G ^ab a bae l{u ^ Uah)t ■ Integrating this 
equation, we find 

a„(t) = a^to) + J2 G U f dt' a ba (t')e^-» (10) 

ab Jt ° 

We consider the situation where the field is initially in 
the vacuum state. We assume that the evolution time 
t — to and the characteristic atomic lifetime r are large as 
compared to the characteristic optical period T. Under 
these conditions, since the continuum of the field modes 
is broadband, the Markov approximation <7f, (i') = a ba {t) 
can be applied to describe the back action of the second 
term in Eq. 1|10|) on the atom [23]. Under the condition 
t — to ^> T, we calculate the integral with respect to t' 
in the limit t — to —* oo. We set aside the imaginary 
part of the integral, which describes the frequency shift. 
Such a shift is usually small. We can effectively account 
it by incorporating into the atomic frequency and the 
surface-atom potential. With the above approximations, 
we obtain 

a^t) = a^to) + 7T Y G*^ ab a ba (t)S(u} - uj ab ). (11) 

ab 

Inserting Eq. (jl 1|> into Eq. © yields the Heisenberg- 
Langevin equation 

d = \ E {laa-bvK'V^bae-^-^'^ 
aa'bb' 

+ l W ^[0,^«-]e iKi - u '' l ' ]i )+e O - (12) 

Here. 

laa'bb' = 27T ^ G* uab G ^ b > S(uj - UJ ab ) (13) 



III. BASIC EQUATIONS FOR SPONTANEOUS 
RADIATIVE DECAY OF THE ATOM 

We use the mode expansion approach and the 
Weisskopf-Wigner formalism [2£j to derive the micro- 
scopic dynamical equations for spontaneous radiative de- 
cay of the atom in the surfaced-induced potential. For 
convenience, we introduce the notation a uu i — |^)(^'| 
and uj uv i = lo v — uj v ' for the transition operator and the 
transition frequency, respectively, of an arbitrary pair of 
levels v and v' . We first study the time evolution of an 
arbitrary atomic operator O. The Heisenberg equation 
for this operator is 

fiab 

+ G; ab al[0.a ba }e^-^ t . (9) 



are the decay coefficients and £e> is the noise operator. 
We emphasize that Eq. (|12|) can be applied to any atomic 
operators. 

We now examine the time evolution of the reduced den- 
sity operator p of the atomic system. We multiply Eq. 
(|12fl with p(0), take the trace of the result, use the re- 
lation Tr[O(i)p(0)] = Tt[O(0)p(t)], transform to arrange 
the operator O(0) at the first position in each operator 
product, and eliminate O(0). Then, we obtain the Liou- 
ville equation 

P = ^ haa'bb'icrbaPVa'b' - 5 bb ><T a > a p)e~ %( - Uab ~' JJa ' b ' )t 

aa'bb' 

+ 7:a> b b'(<?b'a'Paab-5 w pa aa ,)e^~^ t ]. (14) 
For the matrix elements p vu i — {v\p\v'), the above equa- 



5 



tions yields 



Paa' 



Pab 



Pbb> 



a" 
a' 

J $>aaW + T^bK^'^^'/W, (15) 



and, in the case where the dipole moment d eg is parallel 
to the interface, we find 



' aa'bb' 



lf M{ / Re (l+e)Fab(ik ab )F;, b ,^k ab ) 



+ / moiab(zk ab )i a , b ,(tk ab )dA. (17) 



where 7 aa - = J2b~faa>bb = 2TrJ2^ b G ^ab G ^'bS{uJ - w a6 ). 
Equations (|15fl show that the spontaneous decay rate 
for a transition a — > & is given by 7 a (, = jaabb — 
2tt^2^ |G /JQ f ) | 2 5(tJ— cj ab ). Meanwhile, the total decay rate 
for an excited level a is -y a — 7aa = ^J2^b \ G ^ab\ 2 5{uJ - 
LUab)- It is obvious that j a = J2b"fab- 

Equations i|15|) correspond to the incoherent radiative 
decay part of the master equation for an atom interact- 
ing with external fields in the proximity of a surface. A 
set of more general equations is presented in Appendix 
El for the case where the atom is driven by classical co- 
herent plane-wave fields. We note that the terms in Eqs. 
(|15|1 are different from and much more complicated than 
the conventional phenomenological decay terms |2l| . The 
unusual structure of Eqs. I|15fl results from the complex- 
ity of the translational levels of the atom in the surface- 
induced potential. Similar structures are also observed 
in the case of a multilevel alkali atom 4, 221. 



IV. 



SPONTANEOUS RADIATIVE DECAY 
RATES 



In the half-space x > 0, where the atom is restricted 
to, the mode functions of the radiation field are described 
by expressions iJSJ. We insert these expressions into Eq. 
© and then the result into Eq. (|13|1 . We perform the 
summation with respect to the mode index /i. Then, in 
the case where the dipole moment d eg is perpendicular 
to the interface, we obtain 



I aa'bb' 



- lf M{ / Re (l~e)Fab(^ ab )F:, b ,^k ab ) 



r±(OF ab ^k ab )F a , b ,^k ab ) d£ 



+ J T ± (Olab(ik ab )I a , b ,^k ab )dA, (16) 

' 



Here we have introduced the notation 



rx(o = a-a 

m(0 



At- yX-i + f* 



(18) 



and 
T l|(0 



i\ - 1 (nf + \)e + 1 



€(1 



[n\ + l)^ 2 + 1 



The quantity 



7/M 



< n 



(19) 
(20) 



37reo?ic 3 

is the natural linewidth of a two- level atom with the tran- 

leg- 



sition frequency to and the dipole moment d en . The ma- 
trix elements 



F ab (f3) = (Pale^M, 



(21) 



for the transitions between the translational eigenstates 
\ip a ) and \cp b ) have been introduced. In deriving Eqs. 
(|16|l and (|17l) we have changed the integration variable 
K to £ = y|l — k 2 |. Here the parameter k = K/k ab with 
/c a b = uj ab /c has been introduced. 

The functions r± and rn are related to the reflection co- 
efficients r| x = (£-??)/(£ and r 2i = (ni€-v)/(nH + 



rj) as given by r±_ — k rrfi and 1 = 



'21 



Here, 



k = £ 2 an d f? = \/n 2 — 1 + Hence, the terms 

that contain r± and rii in Eqs. (|16|) and (|17|l are the re- 
sults of the interference between the emitted and reflected 
fields. Meanwhile, the functions T± and Tji are related 
to the transmission coefficients ff 2 = 277/(77 + i£) and 
^2 = 271177/(77 + in 2 £) of the evanescent modes as given 
byT ± = (C/2 V W\t p 12 \ 2 andT|| = (£/277)(fe | 2 +^ 2 | 2 ). 
Here, k = \/l + £ 2 and 77 = \/n 2 — 1 — £ 2 . Hence, the 
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terms that contain T± and Tn in Eqs. I|16|l and l|17|) are 

the results of the emission into the evanescent modes. 

The matrix elements F ab ((3) and I a b{P), defined by 
Eqs. 1^21(1 . depend on the overlap between the transla- 
tional wave functions ip a and ip b . The factors |i 7 a b(0)| 2 
and |/ a h(0)| 2 , with the argument (3 = 0, are the same and 
equal to the Franck-Condon factors | (v^o, | | 2 - The ex- 
ponential factors e l/3x and e~@ x in Eqs. I|21|l take into 
account the momentum of the photon emitted or ab- 
sorbed by the atom (see Appendix [Bj. The factor e l/3x 
corresponds to the case where the mode is a propagating 
mode, and the factor e~° x corresponds to the case where 
the mode is an evanescent mode. Thus, F ab {fi) character- 
izes the strength of the translational transition with the 
participation of a propagating- mode photon, and I a b{P) 
corresponds to the case of an evanescent-mode photon. 

We assume that the orientation of the atomic dipole 
moment is completely random. In this case, the aver- 
aged decay parameters are given by j aa 'bb' = {laa'bb' + 
2 7a'a'hf/)A Using Eqs. 1(15)1 and !fT7|l. we find 



laa'bb' = If 



r±(0 + m(0 

2 



F 



\b{£k a b)F a >b' {£,k (lit 



J J±(fl + T ^ Iab{Ckab)Ia , b ,^ kab)d A. (22) 



We note that all the decay coefficients are real quantities. 

When we set a' — a and b' = b in Eq. (|2*2*|) . we obtain 
the following expression for the spontaneous transition 
rate: 



7ab = j f (tOab){ I \ \F ab (£k ab )\ 2 

Re[F^k ab )]}d^ 

I 2 ab(tkab)dA. (23) 



Equation (|23(l shows that the effects of the dielectric on 
the spontaneous transition rate 7 Q b appear through (1) 
the shift of the transition frequency oj ab , (2) the overlap 
\F a b\ 2 between the center-of-mass wave functions, (3) the 
interference between the emitted and reflected fields, and 
(4) the transmission into the evanescent modes. 

The cross decay coefficient j a a' = J2 b laa'bb is the sum 
of laa'bb over the ground-state levels b. The summation 
can be simplified if, for each fixed index a, the overlap- 
ping factors F ab and I ab are substantial only in a small 
region of b. In this case, the transition frequency u ab in 
expression (|22|l can be replaced by an effective frequency 



uj ab that does not depend on b. We use this approxima- 
tion, set b' = b, and sum up the result with respect to b. 
Then, with the help of the relation J2 b I'P^ifbl = 1, we 
obtain 



laa' = {PahxiUab^Va' 



(24) 



Here. 



7x(w) = 7/(cj)<1+ / - — cos(2£kx)d£ 



m+m e -^A (25) 



is the spontaneous emission rate of a two-level atom with 
frequency u>, being at rest in the vicinity of the dielectric 
jl3l Il4j|. In order to find the linewidth 7 a , we set a = a' 
in Eq. J23J. Then, we obtain 



7a = {tp a \l x {uJ ab )\<p a ). 



(26) 



The above equation means that the radiative linewidth 
7 a can be approximately considered as the average of 
the effective- frequency rest-atom decay rate "f x (u> ab ) with 
respect to the center-of-mass wave function ip a (x). 

We note that, in the case where the refractive index n\ 
of the dielectric is not large, the reflection of light from 
the interface and the emission of light into the evanescent 
modes are not strong. In this case, j x (w) with a fixed 
argument w is a slowly varying function of x. Hence, for 
a 7^ a', since |a) and \a!) are orthogonal to each other, 
the quantity 7 aa ' is usually small compared to 7 a and 

la'- 

We emphasize that the first integral in Eq. (|25|l results 
from the interference between the emitted and reflected 
fields. The second integral in Eq. (|25|l results from the 
emission into the evanescent modes. We note that expres- 
sion 1251) for 7a; is in full agreement with the results of the 
linear-response formalism 0,0]. However, the results 
ofRefs. [II El are written in a different form that is 
difficult to identify the origins and physical meanings of 
the contributions. 

If we neglect reflection from the interface and emis- 
sion into the evanescent modes and limit ourselves to 
considering only levels with negligible surfaced-induccd 
frequency shifts, then Eq. (1221) reduces to 

laa'bb' = lofaa'bb' , (27) 

where 70 = 7/(0^0), = w e — u> g , and 

faa'bb' = \ f^i^e^Mi^le-^l^,) 



<lrdx' Smk ? X (p a (x)(p b (x)tp a '(x')(p b '(x'). 



ko(x — x') 



— oc^oo 



(28) 
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Setting b' = b in expressions (EH and (HI and sum- 
ming up the resultant expressions over 6, we find j aa r — 
7o{a\a'). 

V. NUMERICAL RESULTS 

In what follows, we present the results of our numerical 
calculations for the spontaneous radiative decay charac- 
teristics of the atom. As stated in Sec. [n] we use the 
parameters of fused silica, for the dielectric, and the pa- 
rameters of atomic cesium with the Di line, for the atom. 
The refractive index of the medium is n\ = 1.45. The 
wavelength of the D2 line of atomic cesium in free space 
is Ao = 852 nm. Other parameters are given in Sec. [HI 
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FIG. 6: Normalized decay rate 73;(a;)/7/(a;) of a two- level 
atom rested at a point in the vicinity of a semi-infinite di- 
electric medium as a function of the normalized surface-atom 
distance kx. The refractive index of the medium is n\ — 1.45. 

We first consider the case where the center-of-mass 
motion of the atom can be neglected. We plot in Fig. 
the normalized decay rate "f x (uj)/"ff(uj) as a function 
of the normalized surface-atom distance kx. The figure 
shows that "f x (uj)/jf(Lu) varies slowly with kx. We ob- 
serve not only enhancement, 7^/7/ > 1, but also inhibi- 
tion, 72/7/ < 1, of spontaneous emission, depending on 
the position of the atom. Such changes are quantum elec- 
trodynamic effects resulting from the modifications of the 
field modes in the presence of the dielectric 0, Q . The 
maximum value of j x (Lu)/jf(uj) is about 1.6, achieved at 
kx = 0. Such a quantum enhancement is moderate. It 
is not dramatic. The reason is that the refractive index 
of silica rii = 1-45 is not large. We observe small oscil- 
lations in "f x (u))/jf(to) as kx increases. Such oscillations 
are due to the interference between the emitted and re- 
flected fields. 

The deviation of j x (uj)/jf (ui) from unity is caused by 
the interference between the emitted and reflected fields 
and by the emission into the evanescent modes, which are 
expressed by the first and second integrals in Eq. 1)25(1 . 
respectively. We calculate these integrals separately and 
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FIG. 7: Contributions of the first (solid line) and second 
(dashed line) integrals in Eq. (1251 to the normalized decay 
rate y x (u)/jf(u) of Fig. 

plot the results in Fig. [7| The solid curve of the fig- 
ure shows that, for kx < 8, the first integral is negative. 
Thus, in the close vicinity of the interface, the interfer- 
ence between the emitted and reflected fields is destruc- 
tive and hence tends to reduce the spontaneous decay 
rate. The solid curve of the figure also shows that, in the 
region of large kx, the first integral can become positive, 
that is, the interference can be constructive, depending 
on the position of the atom. Meanwhile, the dashed curve 
of the figure shows that the second integral in Eq. I|25|l is 
always positive. Thus, the emission into the evanescent 
modes always enhances the spontaneous decay rate. Such 
an enhancement is due to the existence and confinement 
of the evanescent modes. 

Due to the surface-atom potentials V e and V g , the 
transition frequency of the atom varies in space as given 
by oj — ui x = loq + [V e (x) — V g (x)]/h. Consequently, when 
we take into account the frequency shift, the decay rate of 
the atom is "/ x (lu x ). The spatial dependence of "f x (co x ) is 
determined not only by the change in the mode structure 
but also by the change in the atomic transition frequency 
u> x . We plot the spatial dependence of j x (us x ) in Fig. [S] 
The figure and the inset show that, when x decreases 
from 1 nm to 0.2 nm, the decay rate drops quickly from 
its peak value. The reason is that, in this region of space, 
the surface-induced frequency shift of the transition is 
negative (red shift) and substantial. Such a shift of uj x 
reduces the natural linewidth 7/(1^), see Eq. J20|), and 
hence reduces the decay rate j x (uj x ), see Eq. Q25JI. The 
effects of the transition frequency shift and the quantum 
enhancement on the decay are substantial in the regions 
x < [(C 3e - C 3g )/hu;o} 1/3 = 0.16 nm and x < X = 852 
nm, respectively. The two scales are quite different from 
each other, namely, [{C 3e — Cz g ) /hujo} 1 ^ <C A . In the 
region x > 1 nm, the frequency shift of the transition is 
small and hence the spatial dependence of the decay rate 
is mainly determined by the mode structure of the field. 
Meanwhile, in the region from 1 nm to 0.2 nm, the ef- 
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FIG. 8: Normalized decay rate fxi^x)/^ of a two-level atom 
with a space-dependent frequency, rested at a point in the 
vicinity of a semi-infinite medium. The transition frequency 
uix is shifted from the free-space transition frequency ojo by the 
surface-atom interaction. The parameters of the potentials 
are as in Fig. Q The refractive index of the medium is n\ — 
1.45. The wavelength of the transition of the atom in free 
space is Ao == 852 nm. 



feet of the frequency shift is dominant. We note that, in 
the region x < 0.2 nm, which is not shown in the figure, 
the transition frequency shift may become positive lead- 
ing to an enhancement of the decay rate. In this region, 
the potential slope is steep, the force is large, and con- 
sequently the consideration of the center-of-mass motion 
of the atom must be included. Due to the lack of quanti- 
tative information about the silica-cesium repulsive po- 
tential and due to the need to include the center-of-mass 
motion of the atom, we do not plot the decay rate j x (lo x ) 
in the region x < 0.2 nm. The cutoff value of 0.2 nm is 
chosen because it is close to the position x m = 0.19 nm 
of the minima of the potentials V e and V g . 

In the close vicinity of the surface, we must take into 
account the center-of-mass motion of the atom. To do 
this in a fully quantum treatment, we must consider the 
sets of combined states { | a) } and { | b) } instead of the in- 
ternal states |e) and \g), respectively. We plot in Fig. 
HO the spontaneous transition rate 7 Q0 between shallow 
levels a and b, with large vibrational quantum numbers. 
The figure shows that each upper level a can be substan- 
tially coupled to one or several lower levels b. In general, 
jab varies unevenly with increasing v a or v\>. The values 
of v\, for which j a f, is substantial tend to become larger 
when v a increases. Such features are due to the fact that 
7 Q b substantially depends on the overlap and interference 
between the center-of-mass wave functions ip a and tpb- It 
is clear from Figs. [31and[5]that the frequency shifts of the 
significant transitions between the bound levels a and b 
with large vibrational quantum numbers are small com- 
pared to the free-space optical frequency luq. 

We use the approximate equation (|26|) to calculate 
the linewidths 7 a of shallow bound excited-state levels 
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FIG. 9: Spontaneous transition rates 7 a t for levels with large 
vibrational quantum numbers. The rates are normalized to 
the free-space spontaneous decay rate 70. The parameters 
used are as in Figs. QUI and|H] 
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FIG. 10: Radiative linewidths 7^ of bound excited-state levels 
with large vibrational quantum numbers. The linewidths are 
normalized to the free-space spontaneous decay rate 70. The 
parameters used are as in Figs. and|5] 



a, with large quantum numbers v a , in the range from 
385 to 429, and plot the results in Fig. E3 The fig- 
ure shows that 7 a /7o increases slowly from 1 to 1.53 as 
v a reduces from 429 to 385. The observed enhancement 
of 7 a is the average of the quantum enhancement of the 
rest-atom decay rate ^ x (uj a i) with respect to the center- 
of-mass wave function <p a {x) [see Eq. (|26|) ]. When v a is 
large enough, the effect of the surface-atom interaction 
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on the transition frequency is small, that is, u a i = luq. 
In this case, the magnitude of 7 a is mainly determined 
by the mode structure. When v a is not too large, the 
center-of-mass wave function ip a (x) is confined in a small 
spatial region close to the surface and, hence, the quan- 
tum enhancement of 7 a is observed. The spatial spread 
of ip a (x) increases with increasing v a . This explains why 
"f a reduces with increasing v a in Fig. 1101 When v a is 
very large, the quantum enhancement is negligible and, 
hence, we have 7 a = 70 • 
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FIG. 12: Same as Fig. IHJbut for levels with small vibrational 
quantum numbers. 



FIG. 11: Normalized density y a b f /lo for the rate of sponta- 
neous radiative decay from a bound excited-state level a to 
a free ground-state level bf as a function of the vibrational 
quantum number v a - The parameters used are as in Figs. Q 
El andEl 

The spectra of translational levels of ground and ex- 
cited states include not only bound levels but also free 
levels. Both types of levels can be coupled to each other. 
Since the center-of-mass wave functions of free excited- 



state levels a 



and free ground-state levels bf are nor- 



malized per unit energy, the quantities "fab J ^ 7a/b> and 
7o*6f are the densities of the spontaneous decay rates for 
the transitions a — > bf, af — > b, and af — > bf, respec- 
tively. The knowledge of these quantities is required for 
the study of the dynamical and spectroscopic character- 
istics of moving atoms. For example, in order to calculate 
the excitation spectrum of an atom initially prepared in a 
free ground-state level bf, we need to know the linewidths 
of the transitions a — > bf. In this connection, we plot in 
Fig. 1111 the normalized density j a b, /70 for the rate of 
spontaneous decay from a bound excited-state level a to 
a free ground-state level bf as a function of the vibra- 
tional quantum number v a . 

The effect of the surface-induced shifts of the transi- 
tion frequencies on the transition rates and linewidths be- 
comes substantial when the translational eigenfunctions 
of the atom are close to the surface, that is, when the 
levels are deep. We plot in Figs. El an d El the spon- 
taneous transition rates j a b and the linewidths j a , re- 
spectively, for deep bound levels, with small vibrational 
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FIG. 13: Same as Fig. HOI but for levels with small vibrational 
quantum numbers. 



quantum numbers. Our additional numerical calcula- 
tions (not presented here) show that the transitions from 
bound excited-state levels with small v a to free ground- 
state levels are negligible. Due to this fact, the sum of 
the values of the individual transition rates ~/ a b presented 
in each fixed row of Fig. El almost coincides with the 
corresponding value of the linewidth j a in Fig. 1131 Ac- 
cording to Fig. El m the region of small v a , we have 
la < 70; and ~f a reduces with decreasing v a . In this 
region, the surface-induced transition frequency shift is 
negative and substantial, leading to a decrease in w Q t> 
and, consequently, in jf(u) a t,), 7 a b, and j a . The mini- 
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mum value of the radiative linewidth j a is observed for 
u a = and is about O.2870 This value is substantially 
(about 3.6 times) smaller than the natural linewidth 70 of 
the atom. We note that, when the atomic center-of-mass 
states are close to the surface and the surface tempera- 
ture is high, the decay of the atom due to the interaction 
with the phonons of the solid may become important . 
However, the study of the decay due to phonons is be- 
yond the scope of this paper. 



VI. CONCLUSIONS 

We have studied spontaneous radiative decay of trans- 
lational levels of an atom in the vicinity of a semi- 
infinite dielectric. We have systematically derived the 
microscopic dynamical equations for the spontaneous de- 
cay process. We have calculated analytically and nu- 
merically the spontaneous transition rates between the 
bound excited and bound ground states and the radia- 
tive linewidths of the bound excited states. We have 
shown that the effects of the dielectric on the sponta- 
neous transition rates and radiative linewidths appear 
through (1) the shift of the transition frequency, (2) the 
overlap between the center-of-mass wave functions, (3) 
the interference between the emitted and reflected fields, 
and (4) the transmission into the evanescent modes. Our 
numerical calculations for the silica-cesium interaction 
have demonstrated that the radiative linewidths of the 
bound excited levels with large enough but not too large 
vibrational quantum numbers are moderately enhanced 
by the emission into the evanescent modes and those for 
the deep bound levels are substantially reduced by the 
surface-induced red shift of the transition frequency. We 
believe that the radiative longevity of deep bound states 
can find important applications for atom trapping, quan- 
tum information, and atom and photonic nano-devices. 
We emphasize that we have studied in this paper only 
spontaneous radiative decay. The problem of decay due 
to phonons of the solid is beyond the scope of our paper. 
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be positive and large enough to create a steep slope for 
the potential V(x) in the small interval (x p ,x m ). Such a 
slope leads to a substantial repulsive force on the atom 
in the close vicinity of the surface. We discard the region 
x < x p , assuming that the atom is outside and cannot 
penetrate into this region of space. 

The depth of the potential is defined by D = 
-V(x m ) = -Ae- ax ™ + C3/24. We find the relation 



.4 



3D 



In terms of the parameter x a = y C3/D, we have 

3 



(Al) 



(A2) 



Equations (|A1(I and (|A2(I allow us to determine the pa- 
rameters A and a from the parameters C3, D, and x m . 
Since a > and ax m > 4, the condition 



(A3) 



must be satisfied. 

We use a few available theoretical and experimen- 
tal data to estimate the silica-cesium potential param- 
eters. According to Ref. j2^|, the van der Waals co- 
efficient for the interaction between a ground-state ce- 
sium atom and a perfect metal surface is ^ motal ) = 
4.5 a.u. = 4.39 kHz fim 3 . For a dielectric medium of 
refractive index n, we have an approximate expression 
C 3 = {n 2 -l)/{n 2 + l)C^ motal) 0. For pure fused silica, 
with n = 1.45 (in a broad region around the wavelength 
Ao = 852 nm), we find C 3g = 1.56 kHz ^irn 3 . According 
to Ref. j2^|, the ratio between the van der Waals co- 
efficients for excited- and ground-state cesium atoms is 
C 3e /C 3g = 1.98. This yields C 3e = 3.09 kHz /mi 3 . Re- 
garding the potential depths, we use the experimental 
value D g — 0.66 eV = 159.6 THz, measured as the ad- 
sorption energy of ground-state cesium atoms on fused 
silica j2^| . With the assumption D e /D g — C e /C g , we 
find D e — 316 THz. For both the ground- and excited- 
state potentials, we find x a — 0.21 nm. Hence, the con- 
dition (| A3|) reads 0.13 nm < x m < 0.21 nm. We assume 
that the minimum positions of both potentials are the 
same and are equal to x m — 0.19 nm. Then, Eq. l|AT)l 
yields^ = 1.6 x 10 18 Hz and A e = 3.17 x 10 18 Hz, while 
Eq. (|A2|) gives a g = a e = 53 nm -1 . 



APPENDIX A: PARAMETERS OF THE 
SILICA-CESIUM POTENTIAL 

Consider the potential V(x) = Ae~ ax — C 3 /x 3 . Under 
the condition a 3 Cz/A < 256/3e 4 , this potential has a 
peak at x p = £, p /a and a local minimum at x m = £ m /c«, 
where £ p < 4 and £ m > 4 are the solutions to the equa- 
tion £ 4 e _, » = 3a 3 C 3 /A. The peak value V(x p ) should 



APPENDIX B: DENSITY-MATRIX EQUATIONS 
FOR A COHERENTLY DRIVEN ATOM 

We present the density-matrix equations for a coher- 
ently driven atom. We assume that the atom is driven 
by one or several classical coherent plane-wave fields 
Ej. Here the index I labels the fields. For simplic- 
ity, we consider the case where all the driving fields 
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propagate perpendicularly to the surface of the dielec- 
tric. The expressions for the fields are given by E; = 
[£ ie i(ft*-«it) 6l +c.c.]/2, where £ u e h u u and |$| = wj/c 
are the envelopes, the polarization vectors, the frequen- 
cies, and the wave numbers, respectively. The quanti- 
ties (3i are positive or negative for rightward- or leftward- 
propagating fields, respectively. 

The interaction between the atom and the driving 
fields is described by the Hamiltonian 

Hi = ~ J2( n ^e- iuJit \a)(b\ + H.c), (Bl) 

lab 

where 

n lab = ^F o6 (A) (B2) 

is the Rabi frequency for the action of the field E; on 
the transition between the translational levels \a) and 
\b). Here, di eg — e; ■ (e|d|p) is the projection of the 
atomic dipole moment onto the polarization vector e;, 
and F ab (k) = (ip a \e' lkx \ipb) is the overlapping matrix ele- 
ment for the transition between the center-of-mass states 
ip a and ifb with a momentum transfer of Tik. It is inter- 
esting to note that, when we use the momentum repre- 
sentations ifiaip) and ifbip) for the wave functions ip a (x) 
and ipb(x), respectively, we have 

/oo 
<p* a (p+hk)<p b (p)dp. (B3) 
-oo 

Equation (|B3(I indicates that the overlapping factor 
F a b{k) takes into account the change in the momentum 
of the atom. 



When we apply the Schrodinger equation ihp = [Hi, p] 
to the Hamiltonian I|B1|) and include the radiative decay 
terms, described by Eqs. (JT3J, we obtain the following 
equations for the density matrix p of the atom: 



1,6 

a" 

l,b' l,a' 
a' 

Pw = -|£(<W:6e- i<W ' - ntbPave* 5 ^) 

l,a 

+ +<' a Vby { " hb '- Ua ° ,)t P««'- (B4) 



Here, 5i a b = — u) a + w& is the detuning of the driving- 
field frequency lui from the atomic transition frequency 

^ab- 

We note that Eqs. (|B4|> do not include the interaction 
of the atom with the phonons of the solid. We can take 
into account the decay due to phonons by adding phe- 
nomenological terms |2lj . The coefficients of such terms 
are phonon absorption and emission probabilities and can 
be calculated in second-order perturbation theory j2|| . 
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